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Abstract 

We develop the kinetic theory of Hamiltonian systems with weak long-range inter- 
actions. Starting from the Klimontovich equation and using a quasilinear theory, we 
obtain a general kinetic equation that can be applied to spatially inhomogeneous systems 
and that takes into account memory effects. This equation is valid at order 1/N in a 
proper thermodynamic limit and it coincides with the kinetic equation obtained from the 
BBGKY hierarchy. For N — > +oo, it reduces to the Vlasov equation describing collision- 
less systems. We describe the process of phase mixing and violent relaxation leading to the 
formation of a quasi stationary state (QSS) on the coarse-grained scale. We interprete the 
physical nature of the QSS in relation to Lynden-Belks statistical theory and discuss the 
problem of incomplete relaxation. In the second part of the paper, we consider the relax- 
ation of a test particle in a thermal bath. We derive a Fokker-Planck equation by directly 
calculating the diffusion tensor and the friction force from the Klimontovich equation. 
We give general expressions of these quantities that are valid for possibly spatially inho- 
mogeneous systems with long correlation time. We show that the diffusion and friction 
terms have a very similar structure given by a sort of generalized Kubo formula. We also 
obtain non-markovian kinetic equations that can be relevant when the auto-correlation 
function of the force decreases slowly with time. An interest of our approach is to develop 
a formalism that remains in physical space (instead of Fourier space) and that can deal 
with spatially inhomogeneous systems. 
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1 Introduction 



In the treatment of physical systems in nature, a fundamental distinction must be made be- 
tween systems for which the interaction between particles is short-range or long-range. Systems 
with short-range interactions have been studied for a very long time. They are spatially ho- 
mogeneous, the ordinary thermodynamic limit (N — > +00 with fixed N/V) applies, and the 
statistical ensembles are equivalent for iV — > +00. In recent years, a growing number of phys- 
ical systems with truly long-range interactions have emerged and have been actively studied 
by different groups. An impetus has been given by the conference in Les Houches in 2002 [1] 
which showed the connections and the analogies between different topics: astrophysics, two- 
dimensional hydrodynamics, plasma physics,... Systems with long-range interactions have very 
peculiar properties. They can be spatially inhomogeneous (due to the spontaneous formation 
of coherent structures), the ordinary thermodynamic limit N —>■ +00 with fixed N/V does 
not apply and the statistical ensembles are generically inequivalent. This does not mean that 
statistical mechanics breaks down for these systems, but simply that it must be reformulated 
so as to take into account their peculiarities. Therefore, we must go back to the foundations 
and to the basic principles of statistical mechanics, thermodynamics and kinetic theory. 

In previous papers of this series [2J [3j 0] , we have developed a statistical mechanics and a 
kinetic theory adapted to systems with weak long-range interactions. Our approach heavily 
relies on many important works that have been developed in astrophysics for stellar systems 
[H El E] > in hydrodynamics for two-dimensional vortices [SJ El EH] and in plasma physics [HI [T2] . 
However, the originality of our approach is to remain as general as possible and develop a 
formalism that can be applied to a wide variety of systems with long-range interactions. This 
includes important toy models like the HMF model [13], for example, where explicit analytical 
results can be obtained. However, our approach is more general and aims at showing the 
unity of the subject and the connection between different systems, emphasizing their analogies 
and differences. In this sense, we extend our original approach [TU [10] where we first showed 
the analogy between the statistical mechanics of stellar systems and two-dimensional vortices, 
which are two physical systems of considerable interest. 

In this paper, we consider material particles having inertia and interacting via a weak 
long-range binary potential of interaction u(\r — r'|) in a space of dimension d. In Paper I, 
we have determined the statistical equilibrium states and the static correlation functions in a 
properly defined thermodynamic limit. For attractive potentials, we have shown the existence 
of a critical energy E c (in the microcanonical ensemble) or a critical temperature T c (in the 
canonical ensemble) separating a spatially homogeneous phase from a spatially inhomogeneous 
phase. In Paper II, using an analogy with plasma physics, we have developed a kinetic theory of 
systems with long-range interactions in the spatially homogeneous phase. In Paper III, we have 
studied the growth of correlations from the BBGKY hierachy and the connection to the kinetic 
theory. In the present paper, we further develop the kinetic theory of Hamiltonian systems with 
long-range interactions by starting from the Klimontovitch equation and using a quasilinear 
theory. We derive general kinetic equations that can be applied to spatially inhomogeneous 
systems and that take into account memory effects. These peculiarities are specific to systems 
with unshielded long-range interactions and are novel with respect to the much more studied 
case of spatially homogeneous systems with short-range (or shielded) interactions. However, we 
show that when the system is spatially homogeneous and when memory effects can be neglected, 
we recover the familiar kinetic equations of plasma physics discussed in Paper II. 

1 The case of 2D point vortices, that have no inertia, is special and must be treated separately [15] . 
2 The kinetic equations discussed in [3] have the same form as the Landau and Lenard-Balcscu equations of 
plasma physics except that the Fourier transform of the potential of interaction u p i asma {k) ~ e 2 /fc 2 is replaced by 
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This paper is organized as follows. In Sec. 12.11 we derive a general kinetic equation from the 
Klimontovich equation by using a quasilinear theory. This equation is valid at order 0(1/N) 
in the proper thermodynamic limit iV — > +00 defined in Paper I. It coincides with the kinetic 
equation obtained from the BBGKY hierarchy in Paper III. For iV — > +00, this kinetic equation 
reduces to the Vlasov equation. At order 0(1/N) it takes into account the effect of "collisions" 
(more properly "correlations") between particles due to finite N effects. It describes therefore 
the evolution of the system on a timescale Ntp, where trj is the dynamical time. This general 
kinetic equation applies to systems that can be spatially inhomogeneous and takes into account 
non-markovian effects. However, in order to obtain a closed kinetic equation, we have been 
obliged to neglect some collective effects. This is the main drawback of our approach: a more 
general treatment should take into account both spatial inhomogeneity and collective effects. If 
we restrict ourselves to spatially homogeneous systems and neglect memory terms, we recover 
the Landau equation as a special case. Therefore, the collective effects that we have neglected 
correspond to the effects of polarization taken into account in the Lenard-Balescu equation 
when the system is homogeneous (in plasma physics, they lead to Debye shielding). In Sec. 
12.21 we develop a quasilinear theory of the Vlasov equation [T7J [TSl [TS] in relation with the 
process of violent relaxation [2U1 EE] in the collisionless regime of the dynamics. We derive 
a kinetic equation for the coarse-grained distribution function f(r,v,t) and use this equation 
to describe the problem of incomplete relaxation [22] leading to deviations from the Lynden- 
Bell distribution. We show the analogies and the differences between the quasilinear theory 
of the Vlasov equation used to describe the process of violent collisionless relaxation and the 
quasilinear theory of the Klimontovich equation used to describe the process of slow collisional 
relaxation. In Sec. [31 we consider the relaxation of a test particle in a bath of field particles. 
The relaxation of the test particle is due to the combined effect of a diffusion term and a friction 
term. We derive the diffusion coefficient from the Kubo formula and the friction term from a 
linear response theory based on the Klimontovich equation. Like in the previous sections, the 
originality of our approach is to develop a formalism that can describe spatially inhomogeneous 
systems and that can take into account memory terms. If we consider spatially homogeneous 
systems with short memory, we recover the results obtained in Paper II. However, spatial 
inhomogeneity and memory effects can be important in systems with long-range interactions. 
Therefore, in Sec. H] we derive non-markovian kinetic equations that generalize the standard 
Fokker-Planck equations. We consider explicit applications to self-gravitating systems and to 
the HMF model. 



2 Kinetic equations from a quasilinear theory 

In this section, we obtain a general kinetic equation ([TBI describing the collisional evolution of 
a Hamiltonian system of particles with weak long-range interactions. This equation, derived 
from a quasilinear theory of the Klimontovich equation, is valid at order 0(1/N) in the proper 
thermodynamics limit N — > +00 defined in Paper I. Then, we discuss the analogies and the 
differences with the quasilinear theory of the Vlasov equation developed in [T7J [T5J |TS] to 

a more general potential ii{k) that can take negative values in the case of attractive interactions. This change of 
sign is crucial for the stability of the homogeneous phase and is responsible for instabilities (similar to the Jeans 
instability) for T < T c or E < E c . In the case where the homogeneous phase is stable (for T > T c or E > E c ), 
the potential ii(k) enters explicitly in the Lenard-Balescu equation (11-49) through the dielectric function and 
the results can be different from those obtained with the Coulombian potential u p i asma {k). However, when 
collective effects are ignored, we get the Landau equation (11-40) where the potential of interaction u(k) appears 
only in a multiplicative constant (11-43) controlling the timescale of the relaxation. Note finally that the results 
of the kinetic theory depend on the dimension of space d [16] . 
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describe the process of violent relaxation [2"Ut |2"T] in the collisionless regime. 



2.1 The slow collisional relaxation 

The exact distribution function (DF) of a system of particles in interaction is a sum of Dirac 
functions 

f d (r, v,t) = ^TmcKr - r,(t))5(v - v,,(t)), (1) 

i 

satisfying the Klimontovich equation 

^ + v^-V$^ = 0, (2) 
at or ay 

where $d(r, t) = J u(\r — r'\)fd(r', v', t)dr'dv' is the exact potential created by fd- The Klimon- 
tovich equation ([2]) should not be confused with the Vlasov equation (fl4|) which has the same 
mathematical structure but which applies to the smooth distribution function /. The Vlasov 
equation is valid during the collisionless regime (see Sec. 12.21) while the Klimontovich equation 
is exact and contains the same information as the Hamiltonian equations (1-1). We now de- 
compose the exact distribution function in the form fd = f + $f where / = (fd) is the smooth 
distribution function and 5f the fluctuation around it. Substituting this decomposition in Eq. 
(j2J) and locally averaging over the fluctuations, we get 

where L = vj^ — is an advection operator in phase space constructed with the smooth 

field. Subtracting Eq. ([3]) from Eq. (J2J) and neglecting non linear terms in the fluctuations H, 
we obtain the following equation for the evolution of the fluctuations 

f + «/ = V»£. (4, 

Equations ([3]) and (T4J) form the basis of the quasilinear theory. For spatially homogeneous 
systems, they can be solved with the aid of Laplace-Fourier transforms and they yield the 
Lenard-Balescu equation (see, e.g., [12] and Appendix B of Paper II). In the present work, we 
shall proceed differently so as to treat the case of systems that are not necessarily spatially 
homogeneous and not necessarily markovian. Our method avoids the use of Laplace-Fourier 
transforms and remains in physical space. This yields expressions with a clear interpretation 
which enlightens the basic physics. The drawback of our approach, however, is that it neglects 
collective effects. The formal solution of Eq. (jlj) is 

8f(t) = G(t, 0)5/(0) + jT drG{t, t - r) V<5$(t - r)^ (t - r), (5) 



3 As shown in Papers I and III, the proper thermodynamic limit corresponds to N — > +oo in such a way 
that the coupling constant u* ~ 1/N while the individual mass m ~ 1, the temperature (3 ~ 1, the energy per 
particle E/N ~ 1 and the volume V ~ 1 are fixed. This implies that |r| ~ 1, |v| ~ 1. We also have f/N ~ 1 
and Sf/N ~ 1/vN so that $ ~ u*f ~ 1 and 5$ ~ u*5f ~ 1/yN. With these scalings, we see that the terms 
that we have kept in Eq. (j4|) are of order Sf ~ V^/V and /<5$ ~ y/~N while the nonlinear terms that we have 
neglected are of order SfSQ ~ 1 <C \/~N. We also note that the l.h.s. of Eq. is of order / ~ N while the 
r.h.s. of Eq. ([3]) is of order 5f6$ ~ 1. It would have been more relevant to work in terms of the normalized 
distribution function F = f/N. Then Eq. ([3]) can be rewritten dtF + LF = (1/N)C(F) where the advective 
term is of order O(l) and the collision term is of order 1/N. Therefore, this equation describes the evolution of 
the system on a timescale ~ Ntn- For N — > +oo, it reduces to the Vlasov equation dtF + LF = 0. 
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where G is the Green function associated with the advection operator L and we have noted 
f(t) = /(r, v, £) and 5$(t) = 5$(r, £) for brevity. On the other hand, the perturbation of the 
potential is related to the perturbation of the distribution function through 

- V5$(t) = - / F(l -> 0)<J/i(0^i, (6) 

where refers to the position r and we have noted 5fi(t) = Sf(ri, v 1; t). Therefore, considering 
Eqs. (JSj) and (jSJ), we see that the fluctuation of the field V53>(t) is given by an iterative process: 
V<5$(£) depends on Sf±(t) which itself depends on V<5$i(t — r) etc. We shall solve this problem 
perturbatively in the thermodynamic limit N — ► +oo. To leading order, we get 



(v**^\ = -I A J d^F^l -> 0)G 1 (t,0)G(t,0)(5/ 1 (0)5/(0)) 

+ ri^/ dT / dx i dx 2 F "( 1 ^°) Gr i(*>*- r ) G! (*'*- r ) 
x |f"(2 - 0)<tfA(t - r)*/ a (t - r)} J^(t - r) 

+F»(2 - - r)5f 2 (t - r)>||(t - r) 



(7) 



Now, the fluctuation is exactly defined by 

8f(r, v,t) = ]T m5(r - r,(t))5(v - v,(*)) - /(r, v, t). (8) 

i 

Therefore, we obtain 

(SfiSf 2 ) = (^m 2 5(xi -Xi)5(x 2 -Xj)) + (^m 2 5(xi - Xj)5(x 2 - X;)) 

-<J>*(*i - Xi )/ 2 ) - (]Tm5(x 2 - Xj )/i) + /i/ 2 . (9) 

* j 

To evaluate the correlation function, we average with respect to the smooth distribution 
fi/(Nm) or fifj/(Nm) 2 . This operation leads to 

(5fi5f 2 ) = ^j^hf* + m/*5(xi - x 2 ) - hh - /2/1 + A/a, (10) 



so that, finally, 



(SfiSf 2 ) = mA5( Xl - x 2 ) - i/i/a. (11) 



Substituting this result in Eq. ©, we find that 

(w^) = o)> A + I A ^ dT | dxlF , (1 _ 0)G(t| t _ T) 

x {^(1 - 0)A(t - r) A (t _ r) + ^(o _> l)f(t - r)||(i - r)}, (12) 
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where we have regrouped the two Greenians G and G\ in a single notation for brevity. Finally, 
replacing this expression in Eq. we obtain the kinetic equation 



x {^ (1 ^ 0)/l ^ + ? (0 ^ 1)/ ii} 1 _; (13) 

This is identical to the general kinetic equation (33) of Paper III obtained from the BBGKY 
hierarchy (or from the projection operator formalism |23j). We note that the term of order 1/N 
in the l.h.s. comes from the first term in Eq. (jSj). It corresponds to the mere advection of the 
fluctuations by the smooth field in Eq. (j3J), i.e. ignoring the coupling between the fluctuations 
of the field and the smooth distribution function (r.h.s. of Eq. which gives rise to the 
collision term. 



2.2 The violent collisionless relaxation 

To leading order in N — > +00, the smooth distribution function /(r, v, t) is solution of the 
Vlasov equation 

at or ov 

where $(r, t) = J u(\r — r'\)f(r',v',t)dr'dv' is the smooth potential created by /. The Vlasov 
equation describes the collisionless evolution of the system due to mean field effects only, before 
the cumulative nature of the collisions becomes manifest on a timescale t co u ~ Ntjj or larger. 
Starting from an initial condition which is dynamically unstable, the Vlasov equation coupled 
to a long-range potential of interaction develops an intricate filamentation in phase space at 
smaller and smaller scales. In this sense, the fine-grained distribution function f(r,v,t) never 
achieves equilibrium. However, if we locally average over the filaments, the resulting "coarse- 
grained" distribution function /(r, v,t) will achieve a steady state on a timescale ~ tp. Since 
the Vlasov equation is only valid in the collisionless regime t <C t co u, this corresponds to a 
quasi- stationary state (QSS) that will slowly evolve under the effect of collisions on a timescale 
~ Nto or larger. We can try to predict this QSS in terms of a statistical mechanics of the Vlasov 
equation, using the approach of Lynden-Bell [20] developed for collisionless stellar systems (see 
also [21] ). In the case where the fine-grained distribution function /(r, v, t) takes only two 
values and r] , the statistical equilibrium state maximizes the Lynden-Bell entropy 

Sl.b. = ~ /I- In - + U - - J In ( 1 - ^ J \drdv, (15) 



at fixed mass and energy. This leads to the coarse-grained distribution function 

7= — 2 • (16) 

1 4. e /3%(^+*)-M 

Note that the mixing entropy (IT5|) is formally similar to the Fermi-Dirac entropy and the 
equilibrium distribution ffTB]) is formally similar to the Fermi-Dirac distribution. An effective 
"exclusion principle" , similar to the Pauli principle in quantum mechanics, arises in the theory 
of violent relaxation because the different phase levels cannot overlap. We stress that the 
Lynden-Bell theory is based on an assumption of ergodicity. Indeed, it implicity assumes that 
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the phase elements mix efficiently during the dynamics so that the QSS is the most mixed state 
compatible with the integral constraints of the Vlasov equation. This may not always be the 
case as discussed in the sequel. 

We can try to determine the dynamical equation satisfied by the coarse-grained distribution 
function /(r, v, t) by developing a quasilinear theory of the Vlasov equation. We decompose the 
distribution function in the form / = / + / where / is the coarse-grained distribution function 
and / -C / a fluctuation around it. Substituting this decomposition in Eq. ffT4l) and taking the 
local average, we get 

| +L7= |_V?7, (17) 
where L = vj^ — is an advection operator in phase space constructed with the smooth 



field. Subtracting Eq. (1171) from Eq. (TBI) and neglecting nonlinear terms in the fluctuations, 
we obtain an equation for the perturbation 

| + i/=V*|. (18) 

Equations (ITTj) and ffT8l) are formally similar to Eqs. (j3J) and (TjJ of the previous section but with 
a completely different interpretation. In Sec. 12.14 the subdynamics was played by fa (a sum 
of (^-functions) and the macrodynamics by / (a smooth field). The smooth field averages over 
the positions of the 5-functions that strongly fluctuate. In the phase of violent relaxation, the 
"smooth" field / develops itself a finely striated structure and strongly fluctuates. Therefore, it 
is not smooth at a higher scale of resolution and a second smoothing procedure (coarse-graining) 
must be introduced. In that case, the subdynamics is played by / and the macrodynamics by 
/. The coarse-grained field averages over the positions of the filaments. 

The coupled equations (fTTj) and (|18p can be solved by an iterative procedure similar to that 
developed in Sec. 12.11 and we finally obtain 

df - d /"* f 

lk + L f = 7T-u dT dr 1 dv 1 dr 2 dv 2 —(1^0)G 1 (t,t-r)G(t,t-r) 
at ovv Jo J m 

x ( — (2 -> 0)/(n, v 1; t - r)/(r 2 , v 2 , t - r)^(r, v, t - r) 
+—(2 - l)/(r, v, t - r)/(r 2) v 2 , t - r) |4(ri, v x , * - r) \ . (19) 

To close the system, it remains for one to evaluate the correlation function /(r, v, t)f (ri, Vi,t). 
We shall assume that the mixing in phase space is sufficiently efficient so that the scale of the 
kinematic correlations is small with respect to the coarse-graining mesh size. In that case, 



/(r, v, t)/(n, vi, t) = 44S(t - n)5(v - Vl )/ 2 (r, v, t), (20) 
where e r and e v are the resolution scales in position and velocity respectively. Now, 



~P = (/ - f? = P~ f- (21) 

We shall assume, for simplicity, that the initial condition in phase space consists of patches 
where the distribution function takes a unique value f — r) surrounded by vacuum (/ = 0). In 
this two-levels approximation f 2 = i]q x / = r] f and, therefore, 

/(r.v.tj/Cn.vi,*) = e d r e d J(r - n)«5(v - vi)7fab " 7)- (22) 
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Substituting this expression in Eq. (TiTJj) and carrying out the integrations on r 2 and v 2 , we 
obtain 

(97 — d /"* f F 11 F v 

-± + Lf = e d r e d — / dr dr 1 dv 1 — (1^0)G{t } t-T)—(1^0) 

ot ovv J J m m 

x{fi(Vo -h)§p -/^-7)^} t _; ( 23 ) 

where / = f(r,v,t — r) and /i = /(ri,Vi,t — r). This equation is expected to describe the 
late quiescent stages of the relaxation process when the fluctuations have weaken so that the 
quasilinear approximation can be implemented. It does not describe the early very chaotic 
process of violent relaxation driven by the strong fluctuations of the potential. The quasilinear 
theory of the Vlasov equation is therefore a theory of "quiescent" collisionless relaxation. 

Equation (1231) is similar, in structure, to Eq. ([131) describing the collisional evolution of 
the system with, nevertheless, three important differences: (i) the fluctuating force T{\ — > 0) 
is replaced by the direct force F(l — > 0) because the fluctuations are taken into account 
differently, (ii) The distribution function / in the collisional term of Eq. (Tl3l) is replaced by 
the product /(r/o — /) in Eq. (1231) . This nonlinear term arises from the effective "exclusion 
principle" , discovered by Lynden-Bell, accounting for the non-overlapping of phase levels in the 
collisionless regime. This is consistent with the Fermi-Dirac-like entropy ( fl5|) and Fermi-Dirac- 
like distribution ( fl6l) at statistical equilibrium (iii) Considering the dilute limit / <C i] to fix 
the ideas, we see that the equations (|23|) and (fT3|) have the same mathematical form differing 
only in the prefactors: the mass m of a particle in Eq. (Tl3|) is replaced by the mass r]oefe d of 
a completely filled macrocell in Eq. (|23|) . This implies that the timescales of collisional and 
collisionless relaxation are in the ratio 

— ~ (24) 

tcoll VO e r e v 

Since rj efe d 3> m, this ratio is in general quite small implying that the collisionless relaxation 
is much more rapid than the collisional relaxation. Typically, t nco ii is of the order of a few 
dynamical times to (its precise value depends on the size of the mesh) while t co u is of order 
~ Nto or larger. The kinetic equation (123]) conserves the mass and, presumably, the energy. 
By contrast, we cannot prove an if -theorem for the Lynden-Bell entropy (TT5"|) . Indeed, the time 
variation of the Lynden-Bell entropy is of the form 

Sl.b. = \e d r e d v [ drdvdndv lT =-i- =- f drQ(t)G(t,t - r)Q(t - r), (25) 

2 J f{Vo - f)f AVo - fi) Jo 



Q(t) 



rn 



1 



0,t) 



fi(Vo-h)i%-f(Vo-f) dl1 



(26) 



and its sign is not necessarily positive. This depends on the importance of memory terms. In 
addition, even if Eq. (12"3"1) conserves energy and increases the Fermi-Dirac entropy monoton- 
ically, this does not necessarily imply that the system will converge towards the Lynden-Bell 
distribution (fl6l) . It has been observed in several experiments and numerical simulations that 
the QSS does not coincide with the statistical equilibrium state predicted by Lynden-Bell. This 
incomplete relaxation [22] is usually explained by a lack of ergodicity and "incomplete mixing" . 
In fact, very few is known concerning kinetic equations of the form of Eq. (1231) and it is not 
clear whether the Lynden-Bell distribution (1161) is a stationary solution of that equation (and 
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if it is the only one). As explained in Paper III for the kinetic equation ([TBI describing the 
collisional relaxation, the relaxation may stop because the current J vanishes due to the absence 
of resonances. This argument may also apply to Eq. (1231) which has a similar structure and 
can be a cause for incomplete relaxation. The system tries to approach the statistical equilib- 
rium state (as indicated by the increase of the entropy) but may be trapped in a QSS that is 
different from the statistical prediction (ITS]) . This QSS is a steady solution of Eq. (12"B"|) . or more 
generally (fl~9l) . which cancels individually the advective term (l.h.s.) and the effective collision 
term (r.h.s.). This determines a subclass of steady states of the Vlasov equation (cancellation 
of the l.h.s.) such that the complicated "turbulent" current J in the r.h.s. vanishes. This offers 
a large class of possible steady state solutions that can explain the deviation between the QSS 
and the Lynden-Bell statistical equilibrium state (1T61) observed, in certain cases, in simulations 
and experiments of violent relaxation. Other causes of incomplete relaxation, due to the rapid 
decay of the fluctuations in space and time (leading to a small value of the current), will be 
described in Sec. 12.31 

2.3 The case of stellar systems 

The case of stellar systems is special and deserves a specific discussion. These systems are 
spatially inhomogeneous but, due to the divergence of the gravitational force F(l — ► 0) when 
i - ! — > r, the integral in the r.h.s. of Eq. (123]) can be evaluated by making a local approximation 
which amounts to replacing /(r 1; Vi,t) by /(r, Vi,t). This approximation is justified by the 
fact that the diffusion coefficient diverges logarithmically when r x — > r (see below). We shall 
also make a markovian approximation /(r, v 1; t — r) ~ /(r, v 1; t), f(r, v, t — r) ~ /(r, v, t) and 
extend the time integration to +oo. Then, Eq. ( 1231) becomes 

dl — d f +OD f F M F v 

Tt +Ll= *5ft? I dT j m (1 - M)-(l ^ 0, ( - ,) 

*{7i(»-7i)J£-7to-7)fi}, (27) 

where now / = f(r, v,t) and j\ = f(r,v\,t). Making a linear trajectory approximation 
Vj(t — r) = Vi(t) and Ti(t — r) = — v$r, we can perform the integrations on rj and r like in 
Appendix A of Paper II. This yields the generalized Landau equation 

^ + L7 = vr(2 7 r) 3 e ^A J d^dWk^u{kf^ • ™){fl(Vo ~ 7i)^ - 7fa> - 7)||}- 

(28) 

As a result of the local approximation, the effect of the spatial inhomogeneity is only retained 
in the advective term L in the l.h.s. of Eq. (!28j) . The same approximations are made for 
collisional stellar systems leading to Eq. (11-44) of Paper II. Equation (I28p can also be written 
in the form 

where In A = J + °° dk/k is the Coulombian factor. It exibits a logarithmic divergence at 
small and large scales and it must be regularized by introducing some cut-offs, writing In A = 
ln(L max / L min ). The integral at large scales must be cut-off at the system size L max ~ R (or 
Jeans length) which plays the role of the Debye length in the present context (see Paper III). 
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For collisional stellar systems, the integral at small scales must be cut-off at the Landau length 
Lmin ~ Gm/vf yp with Vf yp ~ GM/R corresponding to a deflexion at 90° of the particles' tra- 
jectory. This yields a Coulomb factor In A ~ IniV. In the present context, the integral at small 
scales must be cut-off at the resolution length e r . Therefore, In A = ln(R/e r ). This implies that 
the timescale of collisional relaxation and the timescale of violent relaxation are in the ratio 

t ncoU m ln(R/e r ) , . 

1 3~3 1 AT ■ ^ 

It is easy to check [25] that Eq. (1291) conserves the mass and the energy, that it monotonically 
increases the Lynden-Bell entropy f[T5"l) (if-theorem) and that its only stationary solution is 
the Lynden-Bell distribution ffTB]) . Therefore, the kinetic equation (|2"9~|) tends to reach the 
Lynden-Bell distribution ([16]) . However, there are several reasons why it cannot attain it: (i) 
Evaporation: for self-gravitating systems, it is well-known that the Lynden-Bell distribution 
( fl6|) coupled to the Poisson equation has infinite mass so that there is no physical distribution 
of the form (I16p in an infinite domain. The system can increase the Lynden-Bell entropy 
indefinitely by evaporating. Therefore, the generalized Vlasov-Landau equation (T2"9l) has no 
steady state with finite mass and the distribution function tends to spreads indefinitely, (ii) 
Incomplete relaxation in space: The turbulent current J in Eq. (1291) . or more generally in 

Eq. (|T9|) . is driven by the fluctuations fi = f 2 of the distribution function (generating the 
fluctuations 5$ of the potential). In the "mixing region" of phase space where the fluctuations 
are strong, the DF tends to reach the Lynden-Bell distribution (fl6j) . As we depart from the 
"mixing region", the fluctuations decay (f 2 — > 0) and the mixing is less and less efficient 
|| J || — > 0. In these regions, the system takes a long time to reach the Lynden-Bell distribution 
( TIT)]) and, in practice, cannot attain it in the time available (see (iii)). In the two levels case, 
we have / 2 = f(r] — f). Therefore, the phase space regions where / — > or / — ► r] do not 
mix well (the diffusion current J is weak) and the observed DF can be sensibly different from 
the Lynden-Bell distribution in these regions of phase space. This concerns essentially the core 
(/ — * Vo) an d the tail (/ — > 0) of the distribution, (iii) Incomplete relaxation in time: during 
violent relaxation, the system tends to approach the statistical equilibrium state ( |T6j) . However, 
as it approaches equilibrium, the fluctuations of the gravitational field, which are the engine of 
the evolution, become less and less effective to drive the relaxation. This is because the scale 
of the fluctuations becomes smaller and smaller as time goes on. This effect can be taken into 
account in the kinetic theory by considering that the correlation lengths e r (t) and e v (t) decrease 
with time so that, in the kinetic equation (I2"9l) . the prefactor e r (t)e v (t) — > for t — > +oo. As 
a result, the "turbulent" current J in Eq. (129]) can vanish before the system has reached the 
statistical equilibrium state (|T6|) . In that case, the system can be trapped in a QSS that is a 
steady solution of the Vlasov equation different from the statistical prediction ([TB|) . Similar 
arguments have been given in [21] on the basis of a more phenomenological kinetic theory of 
violent relaxation. On longer timescale, the encounters must be taken into account. Then 
the system is described by a collisional kinetic Vlasov-Landau equation of the form (111-41). 
This equation conserves the mass, the energy (kinetic + potential) and monotonically increases 
the Boltzmann entropy. The mean field Maxwell-Boltzmann distribution (1-24) is the only 
stationary solution of this equation so that the system tends to reach this distribution on a 
timescale (AT/ In N)to- In practice, however, the convergence to the Boltzmann distribution is 
hampered by the escape of stars and by the gravothermal catastrophe [SI El E]- 
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2.4 Physical interpretation of the QSS 

Based on the preceding kinetic theory, we propose the following interpretation [10J of the QSS 
observed in Hamiltonian systems with long-range interactions: 

1. The QSS results from a process of phase mixing and violent relaxation. This is a purely 
collisionless process driven by the fluctuations of the mean-field potential. It takes place on a 
timescale of a few dynamical times where the Vlasov equation is valid. The QSS is a nonlinearly 
dynamically stable stationary solution of the Vlasov equation on the coarse-grained scale, i.e. 
the coarse-grained DF fQss( r > v ) * s a stable stationary solution of the Vlasov equation. Since 
the Vlasov equation admits an infinite number of stationary solutions, it is not easy to predict 
the one which will be dynamically selected by the process of violent relaxation. 

2. In principle, the distribution /os,g(r, v) of the QSS can be predicted from the statistical 
theory of the Vlasov equation developed by Lynden-Bell |20j. The distribution f L B (r,\) 
depends on the details of the initial condition (in addition to the value of the mass and the 
energy) because of the conservation of the Casimir constraints [21]. The coarse-grained DF 
predicted by Lynden-Bell looks like a sort of superstatistics. 

3. In many cases, the prediction of Lynden-Bell works well [261 [27| l28 | 129]. In certain cases, 
the prediction of Lynden-Bell fails because of the complicated problem of incomplete relaxation 
[22\. The system tends to reach the Lynden-Bell distribution (as implied by the increase of 
the Lynden-Bell entropy) but cannot attain it because the fluctuations of the potential (which 
drive the evolution) fade away before the system has reached the most mixed state. Therefore, 
the incompleteness of the violent relaxation is of dynamical origin. In such cases, the QSS can 
take forms that are different from the statistical prediction, i.e. fQ SS ^ f l.b.- Thus, other 
distributions, that are stable stationary solutions of the Vlasov equation, can emerge. For 
example, the Tsallis distributions [30J are particular stationary solutions of the Vlasov equation 
(polytropes) [31] that can sometimes be reached as a result of an incomplete violent relaxation. 
Several examples have been exhibited where the QSS [321 E31 EH ES] or the transient stages 
of the collisional relaxation [37)1 W(\ I3"5] are remarkably well fitted by Tsallis distributions (see 
the detailed discussion of Paper III [1]). This suggests that Tsallis distributions may represent 
"attractors" of the Vlasov equation in case of incomplete relaxation, for some particular initial 
conditions. However, they are not "universal attractors" 0. Indeed, other distributions have 
been observed that differ both from the Lynden-Bell and the Tsallis distributions. This is clear 
for galaxies in astrophysics that are neither isothermal nor polytropic [10]. There are also cases 
where the system does not reach a QSS and develops instead long-lasting oscillations [UJ W2\ . 
It would be interesting to know whether these different possible behaviours are captured by the 
kinetic equation fT23|) . 

4. Since the Lynden-Bell/ Vlasov approach is restricted to the Boltzmann /i-space, that is 
only a projection of the full Gibbs T-space, one could fear that some fundamental properties 
of the latter (fractal structures, etc..) could be lost in that approach. In fact, we believe 
that the Vlasov equation correctly describes the regime where the QSS appears. Therefore, 
in this regime, all the physics of the problem is contained in the Vlasov equation evolving 
in /i-space. However, the Vlasov equations is a very complicated equation (like the Euler 
equations of turbulence for example). In particular, it can exhibit fractal structures and non- 
ergodic behaviours just as the iV-body system does. Therefore, the Vlasov equation is not in 
contradiction with a complex structure of phase space: the striking features that have been 
observed for the iV-body problem such as QSS [Ml E7J [29l [35] , phase-space holes/clumps [4~T| |4~2] . 

4 Tsallis entropies apply when the phase space of a system is fractal or multi-fractal. The fractal properties 
of the process of violent relaxation are not known. For the HMF model, an interesting regime where Tsallis 
thermodynamics seems to apply |38j has been found above a critical magnetization {39] . 
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anomalous diffusion [SB], non-ergodic behaviours etc. should also be observed with the Vlasov 
equation (except if they are due to finite iV-effects which is also a possibility to consider). 



3 Relaxation of a test particle in a bath 

In this section, we study the relaxation of a test particle in a bath of field particles. Specifically, 
we consider a collection of N particles at statistical equilibrium (thermal bath) and introduce a 
new particle in the system. To leading order in N — > +00, the particle is advected by the mean 
flow in phase space. However, due to finite N effects (graininess) , the test particle undergoes 
discrete interactions with the particles of the bath and progressively acquires their distribution. 
We wish to study this stochastic process. The probability density P(r, v, t) of finding the test 
particle in r with velocity v at time t is governed by a Fokker-Planck equation involving a term 
of diffusion and a term of friction. These results are well-known when the system is spatially 
homogeneous and memory effects can be neglected, as in the case of plasma physics. In the 
present work, we shall develop a method that allows to treat spatially inhomogeneous systems 
and that takes into account non-markovian effects. Our approach is also valid if the bath is 
made of an out-of-equilibrium distribution of field particles that evolves slowly so that it can 
be assumed stationary on a timescale Ntn, which is the typical relaxation time of the test 
particle in the bath. This is the case in particular for one dimensional systems for which the 
Lenard-Balescu collision term vanishes at order 0(1/N). Therefore, any stable steady solution 
of the Vlasov equation does not evolve on a timescale Ntr> [HH ESI II]- 



3.1 Diffusion coefficient 

The increment of the velocity of the test particle between t — s and t due to the fluctuations of 
the force is 

Av"= f ^(t')dt'. (31) 

Jt-s 

After standard calculations (see, e.g., Sec. 4.2 of [B]), the second moment of the velocity 
increment can be written 



Av^Av u \ 1 



2s s 



{s + T)(^{t)T u {t-T))dT. (32) 



We shall assume that the correlation function of the force decreases more rapidly than r _1 
(note, parenthetically, that this is not the case for the correlation function of the gravitational 
force which precisely decreases as r _1 |45j). Then, taking the limit s — > +00, we find that the 
diffusion coefficient is given by the Kubo formula 

, .., ,'Av tl Av v \ , /Av»Av u ^ r+ °° 
= ( > = hm 



r+00 

/ {^{t)T u {t - r))dr. (33) 
Jo 



2At I \ 2s 

On the other hand, after straightforward calculations (see, e.g., Sec. 4.1 of [44J), we obtain 

{T»{t)F{t - t)) = N{F"{\ -> 0, t)F{\ -> 0, t - r)> 

= / dridvi^fl -> 0, i)T v {\ -+0,t- r)-i(ri, vO. (34) 
J rn 
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Therefore, combining Eqs. (1531) and (13~4"|) . we get 

f+oo r 

D^ U = drdndv^il ^0,t)T v {l ^0,t-r) — (ri,vi). (35) 
Jo m 

For a spatially homogeneous distribution, the diffusion coefficient reduces to 

/■+OO J" 

= drdr 1 dv 1 F^(1^0,t)F v (1^0,t-r)—(v 1 ). (36) 

If we neglect collective effects, the force (by unit of mass) created by the field particle 1 on the 
test particle can be written (see Paper I): 



F(l -»■ 0,t) = -im J ku(k)e i]l( - r - ri) dk. 



(37) 



At time t — r, we have 



F(l 0, t - t) = -im J kn(A;)e ik(r( *- T) - ri( *- r)) rfk. (38) 

To leading order in iV — > +oo, the particles follow rectilinear trajectories so that i"j(t — r) = 
rj — VjT where = r-j(t) and Vj = Vj(t) denote their position and velocity at time t. Then, we 
get (with x = r — ri and w = v — vi): 

F(l 0, t - t) = -im J kn(A;)e ik(x - wr) dk. (39) 

Substituting this expression in Eq. ( 1361) and carrying the integrations on ri and r, we obtain 
after straightforward calculations 

D" u = n(2n) d m J fc" k v u{kf 5 (k ■ w)f{v 1 )dkdv 1 . (40) 

If we take into account collective effects (see Appendix [B]) , we have to replace Eq. (139]) by 



Then, we get 



F(l -> 0, t - r) = -im I k— e ik(x ~ WT) cik. (41) 

e(k,k-vi) 



D"" = n(2n) d m / k»k u — ^-— 6(k ■ w)f(v 1 )dkdv 1 . (42) 



The calculation of the diffusion coefficient tensor D^ v for different potentials of interaction and 
different dimensions of space is performed in Paper II and in [16J. For one dimensional systems, 
we have the simple result 

^"^wjTiS 1 *" (43) 

where we have used 5(k(v — vi)) = (l/\k\)S(v — Vi) to perform the integration on v\. 
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3.2 Friction coefficient 



In addition to its diffusive motion, a test particle evolving in a bath of field particles undergoes 
a dynamical friction. The friction corresponds to the response of the field particles to the 
perturbation caused by the test particle, as in a polarization process. The test particle modifies 
the distribution of the field particles and the retroaction of this perturbation on the test particle 
creates a friction. The expression of the friction force can be derived from a linear response 
theory starting from the Liouville equation as done in Kandrup [16]. In this section, we show 
that it can also be obtained from the Klimontovich equation. This will make a close connection 
to the quasilinear theory developed in Sec. [2J 

The introduction of a test particle in a bath of field particles modifies the distribution 
function /(r, v, t) of the bath. Since this perturbation is small, it can be described by the 
linearized equation 

W ' + «/ = V»j£ (44) 



dt dv 



whose formal solution is 



Sf(t) = j drG(t, t - r)V5<5>(t - r)^(f - r). 

J 



(45) 



We have used the fact that, initially, Sf(0) = 0. On the other hand, the perturbation of the 
force in r is given by 

- V5$(r, t) = - [ F(l - O^frWdx! + [ T{1 - 0)5(1-! - r P (t))d ri , (46) 
m J J 

where rp(t) denotes the position of the test particle. The second term is the force created by 
the test particle and the first term is the fluctuation of the force due to the perturbed density 
distribution of the field particles. Substituting Eq. ([4"5]) in Eq. ([46]) we obtain 

- V5$(r, t) = i jT dT J rfx i F ( 1 -> °) G &> 1 ~ r )^^ - r )|^(* " r ) 

+ / V(l -> 0)5(n - r P (t))dn. (47) 



This is an integral equation for — V5$(r, t). For a spatially homogeneous system, one can solve 
this equation exactly by using Laplace-Fourier transforms. This is how the dielectric function 
enters in the problem (see Appendix [B]). In order to treat more general systems that are not 
necessarily homogeneous, we shall make an approximation which amounts to neglecting some 
collective effects. We solve Eq. ( 1471) by an iterative process: we first neglect the first term in 
the r.h.s. of Eq. ( [471) keeping only the contribution of the test particle. Then, we substitute 
this value in the first term of the r.h.s of Eq. ([47]) . This operation gives 

- V8<$>(r,t) =-— [ dr f d Xl dr 2 F(l -> 0)Gi(t,t - t)T v (2 -> 1) 
m J J 

x^(t - r)6(r 2 - r P (t - r)) + J T{\ -* 0)5^ - r P (t))dn. (48) 

This quantity represents the fluctuation of the field in r caused by the introduction of the 
test particle in the system and taking into account of the retroaction of the field particles. 
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If we evaluate this expression at the position r P of the test particle and subtract the second 
term (self- interaction), we obtain the friction force felt by the test particle in response to the 
perturbation that it caused. Denoting now by the position of the test particle, we find that 
the friction is given by 



F Poi = ~- J o dr J W^(l-0,i)r(0-l 1 t-r)^(r 1 (t-r) 1 v 1 (t-r)). (49) 
For a thermal bath, where the distribution of the field particles is given by /(ri,vi) = 

Ae -^m(^/2+#(r 1 )) ) we obtain 

F £* = Pf dT / dr^F^l ^ 0,^(0^1,1-^-^-^ f^vj, (50) 

where we have used f(ri(t — r), Vi(t — r)) = /(ri(t), vi(t)) since / is a stationary solution of 
the Vlasov equation. This is equivalent to the result of Kandrup [16] based on the Liouville 
equation but it is obtained here in a simpler manner from the Klimontovich equation. We can 
also obtain this result in a slightly different way. We approximate — V5$(r, i) in Eq. f|44|) by 
the force T[P — > 0) created by the test particle only so that 

«/ . Uf =-F(P^V- (51) 



dt * ' <9v 

This equation can be solved with a Green function yielding 

ft 



Sf(t) = - [ drG(t,t - t)T(P -»• 0,t - r)^-{t - r). 
Jo ^ v 



(52) 



This represents the perturbation of the distribution function of the field particles caused by the 
introduction of a test particle in the system. This perturbation produces in turn a force which 
acts as a friction on the test particle (by retroaction). If we substitute Eq. (152]) in the first part 
of Eq. ([15]) and evaluate this quantity at the position of the test particle, we recover Eq. (T4"9l 
for the friction. 

If we now consider a spatially homogeneous distribution of field particles, the expression of 
the friction force becomes 

Fp "° l = ml dT J ^i^i^'d^O^^a^Ct-r)^^), (53) 

where we have used V\(t — r) = vi(£) to leading order in iV — > +oo. Taking the limit t — > +oo, 
we get 



vol m 



- f +C ° dr [ dr 1 dv 1 F' l (1^0,t)F t '(l^Q,t-T)%(v 1 ). (54) 
■n Jo J ov u 

This is a sort of generalized Kubo relation involving the gradient of the distribution function in 
velocity space instead of the distribution function itself. The nice similarity in the expressions 
of the diffusion coefficient IfSh}) and friction force (54\ ) is worth mentioning. The integrals on 
ri and r can be calculated in the same manner as in Sec. 13.11 and we obtain 

F£ ol = 7r(2vr) d m f d^d^kfk^k^^ ■ w)|4- ( 55 ) 
J & v i 
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In order to take into account collective effects, we can follow the approach of Hubbard [UJ . 
The force (by unit of mass) created in r by the introduction of the test particle is 



F(P -> 0) 



-im 



u{k) 



e(k, k ■ vp 



. e ik(r-r P ) dk) 



(56) 



where the dielectric function takes into account the response of the whole system. The bare 
force due to the test particle alone is 



F(P -> 0) 



-vm 



ku(£0e ik{r - rp) dk. 



(57) 



If we subtract Eq. (1571) from Eq. (1561) . we get the force created in r by the perturbation of 
the distribution function of the field particles caused by the introduction of the test particle. 
Evaluating this force at the position of the test particle, we obtain the friction that it experiences 
as a result of the polarization process 



pol 



-im 



ku{k) 



1 



|_e(k,k-v) 



1 



dk. 



This can also be written 



Fp i = m J ku(k)lm 



1 



e(k,k - v 

Using the identity (B12) of Paper II, we finally obtain 

u(k) 2 



dk. 



r pol 



ir(2n) d m / dv x dk 



:(k,k 



Pfc"5(k ■ w 



(58) 



(59) 



(60) 



If we neglect collective effects and take |e(k, k • v)| 2 = 1, we recover Eq. (j55|) obtained in a 
different manner. Now, the friction force is due not only to the polarization but also to the 
variation of the diffusion coefficient with the velocity of the test particle v. As a result, the 
complete expression of the friction force is 



friction 



r pol ^ 



Q D nu 

dv v 



(61) 



The second term is obtained when we take into account the influence of the fluctuations of the 
force in the trajectory of the test particle, i.e. when we go beyond the rectilinear trajectory 
approximation. As shown by Hubbard [17], this is necessary for the calculation of the friction 
while this is not necessary for the calculation of the diffusion coefficient. From Eqs. (pf2"l) and 
([60]) we get 



friction 



m 



d^dkk^k^^kfh 











5(k 



w 



dv v dvi J |e(k, k ■ v 



(62) 



where we have used an integration by parts in Eq. ( 1601) . When we ignore collective effects, 
expressions (1421) and (1601) for the diffusion coefficient and the friction force can be obtained 
directly from the Hamiltonian equations, by making a systematic expansion of the trajectory 
of the particles in powers of 1/N in the limit iV — > +oo as shown in Appendix IA1 

For a thermal bath, corresponding to the case where the field particles are at statistical 
equilibrium, the distribution function is the Maxwell-Boltzmann distribution 



d/2 



pe 



v 2 



(63) 
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Inserting the identity 

Of 



-/3m/ lVl , (64) 



in Eq. f[6"Ul) . using the 5-function to replace k ■ vj by k ■ v, and comparing with Eq. (|42p . we 
find that 

Koi = -pmD^v v . (65) 

This can be viewed as a generalized Einstein relation. We note that the diffusion coefficient 
and the friction coefficient depend on the velocity of the test particle. We also note that the 
Einstein relation is valid for the friction force F po i due to the polarization, not for the total 
friction force We do not have this subtlety for the ordinary Brownian motion where the 

diffusion coefficient is constant. 

We now consider an arbitrary (steady) distribution of the bath. If we neglect collective 
effects and use Eqs. ( 140!) and (|55|) we obtain after simple manipulations (see Eq. (16) in [IB"]): 

■w = 1 (66) 

Therefore, 

F friction 2Fp G ;. (67) 

We note that the friction force calculated by Kandrup [46] corresponds to the polarization part 
Fp i while Chandrasekhar [48J computes the full friction Ff rict i on . This explains why there is a 
factor 1/2 between their results for equal mass particles (see [46], pp. 446). 
Finally, for ID systems, we have the simple result 

F " =4,vw rs*' <68) 

This expression is valid for an arbitrary (steady) distribution of the bath and it takes into 
account collective effects. Comparing Eq. (1681) with Eq. (143]) . we find that the friction force is 
related to the diffusion coefficient by the relation 

/ n d In f , . 

F pol = D{v)—^. (69) 

This can be viewed as a generalization of the Einstein relation for an out-of-equilibrium distri- 
bution of the bath. 



3.3 The Fokker-Planck equation 

Assuming that the system is spatially homogeneous, the probability density P(v,t) of finding 
the test particle with the velocity v at time t is governed by a Fokker-Planck equation of the 
form 

dP = 1 d 2 f (A^A^) \ _ _d_ f (A^)\ 

dt 2dv»dv»\ At J dv»\ At J' 1 ' 

This Fokker-Planck approach assumes that the stochastic process is markovian (see Sec. [4] for 
generalizations). It also assumes that the higher order moments of the increment of velocity Av 
play a negligible role. This is indeed the case in the N —>■ +oo limit that we consider since they 
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are of order 0(N 2 ) or smaller. At order 0(N x ), we have found that the second (diffusion) 
and first (friction) moments of the velocity increment of the test particle are given by 

(Av^Av v ) ^ (At;") _ dDr | ^ 



with 



2At At aw 1 ' 



= 7r(27r) d m / fcT . ,, Q g(k ■ w)/(v 1 )rfkdv 1 , (72) 

e k, k • v r 



es i£ = / |e( ff )2 v)|2 ^(k • w)||rfv lrf k. (73) 

Note that we have changed the sign of 7/ M with respect to Paper II. The Fokker-Planck equation 
rOl) can be written in the alternative form 

dP d ( dP \ 

1 D^^--Pif). (74) 



dt dv^ \ dv 

The two expressions (1701) and (1741) have their own interest. The expression (1701) where the 
diffusion coefficient is placed after the second derivative d 2 (DP) involves the total friction force 
Ff motion = (At> /1 )/At and the expression (174"|) where the diffusion coefficient is placed between 
the derivatives dDdP isolates the part of the friction rf = F^ ol due to the polarization. This 
alternative form (1741) has therefore a clear physical interpretation. Inserting the expressions 
(172]) and (173]) of the diffusion coefficient and friction term in Eq. (1741) . we obtain 

§ = / ^^Rif^^'^-^)!^-^)^'^' 4 '- < 75 > 

For a thermal bath, using Eqs. ( 1651) . the Fokker-Planck equation ( 1741) can be written 

<9P 9 



<9t dv^ 



dP 

D tM/ (—+ (3mPv u 



(76) 



dP _ d 
dt dv 



where D^ v {v) is given by Eq. (1721) . Since the r.h.s. of Eq. (1751) is of order 0(1/N), the 
distribution of the test particle P(v,t) relaxes to the Maxwellian distribution on a typical 
timescale Ntp (see [TB] for more details). In one dimension, the bath f(v) can be any stable 
stationary solution of the Vlasov equation. Using Eq. (I69p . the Fokker-Planck equation (1741) 
can be written 

where D{y) is given by Eq. (1431) . The distribution of the test particle P(v,t) relaxes to the 
distribution of the bath f(v) on a typical timescale Nt D [16J. 

In Paper II, we have obtained the Fokker-Planck equation (I75p from the Lenard-Balescu 
equation (11-49) by replacing /(v,t) by the distribution of the test particle P(v,t) and f(y\,t) 
by the static distribution of the bath /(vi). This procedure transforms an integrodifferential 
equation (11-49) in a differential equation (I75|) . The expressions (I7T]) - (I73|) of the diffusion and 
friction were then obtained by identifying Eq. (1751) with the Fokker-Planck equation (1701) . In the 
present paper, we have proceeded the other way round by first determining the moments (ITT|) - 
( T73|) . then inserting them in the Fokker-Planck equation (ITOj) . Note that Hubbard [47] derived 
the expressions (!7T!) - (!73l) of the diffusion coefficient and friction force but did not make the 
calculations explicitly until the end. In particular, he did not explicitly wrote down the kinetic 
equation (1751) that is related to the Lenard-Balescu equation (11-49) discovered independently 
at the same period 
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3.4 The Fokker-Planck equation at T = 



In Paper II and in [16J, we have given various expressions of the Fokker-Planck equation (175]) 
for different potentials of interaction and different dimensions of space. However, we have 
not explicitly considered the case T = which presents interesting features. At T = 0, the 
Maxwell-Boltzmann distribution (153"]) reduces to /(vi) = p8(vi). Substituting this expression 
in Eq. (172]) . we find that the diffusion coefficient becomes 

= n(2n) d pm [ fcT^^5(k • v)dk, (78) 
J l e C k . 0)1 

with e(k, 0) = 1 + (2Tr) d u(k)/3pm according to Eq. (11-13). Note that, in this section, we 
consider the case of repulsive potentials with u(k) > so that the homogeneous phase is stable 
even at T = (see Paper I). We now observe that the integral in Eq. ( 175]) is similar to the one 
in Eq. (11-41). Therefore, it can be written 

Z^ = ^(V-^Y (79) 



v \ v 2 



where 



K d = X d pm / k 



-oo 



u(k) 1 2 



o [1 + (2ir) d u(k)f3pm 



dk, (80) 



with A3 = 8tt 5 and A2 = 8tt 3 . For the Coulombian potential, we have (2iT) 3 u(k)(3pm = k 2 D /k 2 
where ko is the Debye wavenumber (see Paper I). Therefore, the collective effects encapsulated 
in the dielectric function in the denominator of Eq. (IHDl regularise the integral for k — > (this 
is a particular case of the Lenard-Balescu equation). On the other hand, noting that D^ v v u = 
according to Eq. (1791 . we find that the friction force (155]) vanishes. Therefore, at T = 0, the 
Fokker-Planck equation (175]) can be written 

dP d f5^v 2 -v»v u dP \ 

— = K d —{ , — 1. (81^ 



dt dv^ \ v 3 dv L 

This equation admits an infinity of stationary solutions. Indeed, since D^ u v v = 0, any distribu- 
tion P = P(v) depending only on the modulus v = |v| of the velocity is a stationary solution 
of Eq. ([ST]) . Therefore, at T = 0, the test particle does not necessarily relax to the distribution 
of the bath /(v) = pS(v). On the other hand, for one dimensional systems, Eq. (1HT]) reduces 
to dP/dt = so that the distribution of the test particle does not evolve in time. 



3.5 More general kinetic equations 

It is instructive to compare the Fokker-Planck equation (175]) with the more general equation 
obtained from the projection operator formalism [23] . When collective effects are ignored, this 
equation can be written 

dP OP dP d T 1 f 

x {^(1 - + ^"(0 - l )£p} p ^ v ' 1 ~ T >^< r - (82) 

It can be obtained from Eq. ( TTSl by replacing f(v,t) by P(v, t) and f(vi,t) by /(vi). This 
is a sort of generalized "Fokker-Planck" equation involving a term of "diffusion" and a term of 
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"friction". However, strictly speaking, Eq. (1521) is not a Fokker-Planck equation because it is 
non-Markovian. We also note that the "diffusion" term appears as a complicated time integral 
of the force correlation function involving P(r,v,t — r). This can be seen as a generalization 
of the Kubo formula (1351) . Similarly the "friction" force is a generalization of the expression 
obtained in Eq. (1491) with a more complicated time integral. If we consider a thermal bath 
where the distribution of the field particles is the Boltzmann distribution, we get 



i + v f + < F >f = jf *■ / - °> g < m - t) 

x<jV(l -> 0) ■ ^- -PmF(0 -> 1) •vi}p(r,v,t-r) — (ri,vi). (83) 



<9v J ' ' m 

If we come back to Eq. (|82|) . make a Markovian approximation and extend the time integration 
to infinity, we get 

dP dP dP d f +OD f 

m + v * + < F> ^ = ft? 1 *■ Z *" iv ' F "< 1 - 0)G(( - ( - t > 

xl^l-O^+^-D^l^.v.oArLV,), (84) 

where we recall that the coordinates appearing after the Greenian must be viewed as explicit 
functions of time rj(t — r) and Vj(t — r) (see Paper III for more details). For a spatially 
homogeneous system, Eq. ( 1521 takes the simplest form 



f = £ jf * / dridVlF " (1 °' ()F " (1 - °- ( - T) - sr) p(v - ' " r) ^ (Vl) - (85) 

where we have used Vj(t — r) = v, for a spatially homogeneous system. We shall come back to 
this non-markovian equation in Sec. HI If we now make a Markovian approximation P(v, t — 
r) ~ -P(v, t) and extend the time integral to infinity, we get 

This is a Fokker-Planck equation which can be put in the form (lT4"j) with a diffusion coefficient 

= dr / dr 1 dv 1 F fi (1^0,t)F v (1^0,t-r)— (vA (87) 
Jo J rn 

and a friction force due to the polarization 

^ = ~j °° dr j dv 1 dv l F»{l^Q,t)F»{Q^l,t-T)^L{v l ). (88) 

These expressions agree with Eqs. (136|) and (1541 obtained directly from the equations of motion. 
After integration on r and r 1; we recover the Fokker-Planck equation (1T5]) with the expressions 
(172]) and (173")) of the diffusion coefficient and friction term (with |e(k, k- v)| 2 = 1 since collective 
effects are neglected here). 
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4 The non-markovian equation 
4.1 General results 

In this section, we study in more detail the non-Markovian equation (I85l) . If the field parti- 
cles are at statistical equilibrium (thermal bath), using the identity (l64l) . the non-markovian 
equation (1851) takes the form 

^ = A jT dr J dndvxF^l - 0, f)F"(l - 0, f - r)^(v0 + /3m<) P(v, t - r). 

(89) 

It can be rewritten 

where we have introduced the notations 

C^(t) = (F»(t)F u (t - r)) = N{F"{1 -> 0, f)F"(l -»• 0, * - r)), (91) 

W^(r) = N(F»(1 -> 0,t)F"(l -> 0,t - r)<). (92) 

These quantities can be calculated by making the linear trajectory approximation. The first 
quantity has already been studied in Paper II. It represents the temporal correlation of the 
force acting on the test particle. It can be written 

C^(t) = (2ir) d m J Ffc^A;) V ik(v ~ Vl)T /(vi)dv x dk. (93) 

Performing the integration on vi, we get 

C^{t) = (2n) 2d m [ k»k u u{k) 2 e~ ikVT f(\sr)dk, (94) 



where / is the Fourier transform of /. For a Maxwellian distribution of the field particles 
(thermal bath), we have 

C^(r) = (2n) d pm J k^k v u(k) 2 e-^ T e^ 2 IWm dk. (95) 

On the other hand, the function W(r) is given by 

W{t) = m(2n) d J (k • v 1 )Pw(A;) 2 e^ ik(v - Vl) V(vi)rfv 1 rfk. (96) 

Performing the integration on v 1; we get 

W"(t) = -im(2Tt) 2d J dkk^u(k) 2 e- lkvr -^-f(kr). (97) 

For a Maxwellian distribution of the field particles, we obtain 

W»{t) = -ipm(27r) d J dkk^u(k) 2 e- ikvT -^-e- k2r2/2Pm , (98) 

so that, finally, 

W(t) = t(27r) d p^ J k»u(k) 2 k 2 e- i]iVT e- k2T2/2 P m dk. (99) 
Let us now apply these general results to some specific systems. 
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4.2 Self-gravitating systems 



For the gravitational interaction, we can easily perform the integrations in Eqs. (|95|) and (1991) 
by introducing a spherical system of coordinates with the z axis in the direction of v. The 
correlation function C^ u (t) is given by Eqs. (11-94), (11-95) and (11-96). On the other hand, 
after some calculations, we find that 



. , AirpmG 2 , . 
W(r) = — G(z)v, 



VT 



(100) 



where x = (/5m/2) 1//2 v and G(x) is the function defined by Eq. (11-75). Comparing this 
expression with Eq. (11-95), we find that 



W(r) = C\\(v,t)v = C^(v,t)v u . 
Therefore, for the gravitational interaction, we have the equality 

0, t)F u (l 0, t - r)<) = -> 0, t)F"(l -> 0, t - t))v v 



(101) 



(102) 



We stress, however, that this equality is not true for any potential. Using the relation fllOip . 
we can rewrite the non-Markovian equation fl90|) in the form 



dP 
~dt 



d 



_d_ 

dv v 



+ (3mv v P(v,t-r). 



(103) 



For a spherically symmetric system, the distribution P(v, t) depends only on the modulus 
|v| = v of the velocity and we obtain 



dP _ 1 d 

dt v 2 dv 

where (see Paper II): 



r 2 I (hC\ (t. r)[ ^- + (3mv )P(r.t- r) 



(104) 



Cu(r,v) 



AnprnG^ 



VT 



G{x) 



(105) 



If we make a markovian approximation P(v,t — r) ~ P(v,t) and extend the time integration 
to +oo, we recover the Kramers- Chandrasekhar equation 



dP 
~dt 



v 2 dv 



fdP 

v 2 D ][ \{v)[— + !3mPv 



(106) 



with 



DAv) 



n i \a AirpmG 2 

C\\{T,V)dT = ^{ x ) 



~°° dT 

T 



(107) 



This expression exhibits the well-known logarithmic divergence of the diffusion coefficient which 
appears here in the time integration (see a discussion of this issue in Paper II). The divergence 
for t — > is related to the linear trajectory approximation and could be cured by a more 
accurate treatment of binary collisions. The divergence for t — > +oo is more serious. One 
usually introduces a cut-off but this procedure is relatively ad hoc. Alternatively, one could 
consider the non-Markovian equation f ll04p -( fT05l) which is well-posed for any time t. 
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4.3 The HMF model 



As discussed previously, non-markovian effects can be important for self-gravitating systems 
because the temporal correlation function of the force decreases algebraically, like t~ l . For 
neutral plasmas, the situation is different because of Debye shielding. In that case, collective 
effects cannot be ignored in the computation of the force auto-correlation function and they 
are taken into account through the dielectric function in Eq. (11-98). Then, the temporal 
correlation function is given by Eqs. (11-100), (11-109), (11-110) and (11-21) of Paper II, and it 
decreases exponentially rapidly. In that case, the Markovian approximation is valid. Collective 
effects are also important for the HMF model [5TI 152"], When collective effects are ignored, 
it is found that the temporal decay of the correlation function of the force is gaussian, see Eq. 
(11-113). By contrast, when collective effects are properly accounted for, it is found that the 
correlation function decreases exponentially rapidly, see Eq. (11-111). Furthermore, the decay 
rate tends to zero for T — > T c implying a slow decay of the correlations. This may unveil a 
failure of the markovian approximation close to the critical temperature. For that reason, it 
may be useful to derive non-markovian kinetic equations which take into account collective 
effects. 

Collective effects can be taken into account in the non-markovian equation (1901 by making 
the substitution 

fl « 2 - Hif^TF- (108) 

in the expressions (1931) and (|96|) . The correlation function of the force is now given by 

CHr) = ( 27 d ™ f ,. U ?^ ,, e' tk(v ' Vl) V(vi)^v^k. (109) 

J |e(k,k- vi)| 2 

This can be written 

C^{t) = {2Tx) d m J ^A;^(A;) 2 e- ikvT g(k,r)rfk, (110) 

where the function Q(h, r) is defined by Eq. (11-99). On the other hand, the function W(r) is 
given by 

W{t) = m(2n) d [ (k • vi)F 1 - ?r ^2-^e- ik(v - Vl)T /(vi)dv 1 dk. (Ill) 



|e(k,k- Vl )| 2 



This can be rewritten 



WHt) = -im(2ir) d J d\k^u(k) 2 e- ikvT -^-Q(k,r). (112) 

For a Maxwellian distribution of the field particles, the large time asymptotics of Q(h,r) is 
given by (11-109). Using Eqs. f 1 11 j) and (I112p . we can then obtain the large time asymptotic 
of C^{t) and W{t) for r -> +oo. 

Let us now specifically consider the HMF model where the potential of interaction is trun- 
cated to one Fourier mode. For this system, using Eqs. fill 01) and (j!12p . the non-markovian 
equation (1901 can be written 



2 a rt 



OP k 2 d . , 

dt At: dv J 



d 

Q(t) cos{vt)- f3Q (r) sin(wr) 

ov 



P{v,t-r), (113) 
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where Q(r) behaves like 



2N l/2 



Q{r) ~ P [-5) ^^e" 7T , (114) 

W 71^(7^72)1 

for r — > +00. The damping rate 7 and the function -F(x) are defined in Paper II. As discussed 
above, the exponential relaxation time 7 _1 (T) diverges for T — > T c so that the Markovian 
approximation may not be correct close to the critical point. This may be an interesting 
situation to analyze in deeper detail with the non- markovian equation (I113p . 
If we neglect collective effects, we find that 

Q(k, r) = (2n) d f(kt) = pe - k ^/ 2 ^, (115) 

where the second equality is valid for a Maxwellian distribution of the field particles. For the 
HMF model, we obtain 

Q{r) = pe - 2 W. (116) 

This yields a Gaussian decay of the correlations instead of an exponential decay in Eq. (11141) 
when collective effects are accounted for [52]. With the expression ( 1114ft for Q(t), the non- 
markovian equation (190]) becomes 



dP _ pk 2 d 
dt An dv 







, \ 9 . . . 

COSiVT)— — h T smlVT) 
ov 



P{v,t-r). (117) 



We note that for the HMF model, the equality (11021) is not satisfied. 

The previous equations assume that the distribution of the bath is maxwellian. More 
generally, if we come back to the non-Markovian equation Eq. ( 185]) and perform the integration 
on r l5 we obtain 



f = rn{2 ^ [ dT I dv ^ v W^? cos(k ' wr) - £r) p(v ' ( " r)/(Vl) ^ 



(118) 

For the HMF model, this equation reduces to 

dP k 2 d f* f cos(wt) f d d 



at 4^y rf V d "'Rif#U-^J p ("' ( -^W- (119) 

This equation is valid for any steady distribution of the field particles, not only for the statistical 
equilibrium state (thermal bath). 



5 Conclusion 

In this paper, starting from the Klimontovich equation and using a quasilinear theory, we have 
developed a kinetic theory for systems with weak long-range interactions. We have obtained 
general equations that take into account spatial inhomogeneity and memory effects. These 
peculiarities are specific to systems with unshielded long-range interactions. However, in order 
to obtain closed kinetic equations we have been obliged to neglect some collective effects. 
These collective effects can be taken into account for spatially homogeneous systems with short 
memory time (with respect to the slow collisional relaxation time). In that case, we recover 
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the Lenard-Balescu equation of plasma physics with slight modifications (see Paper II and 
footnote 2). It would be valuable to develop a formalism that takes into account both spatial 
inhomogeneity and collective effects. This has been partly done in Paper III, where we have 
obtained two coupled Eqs. (III-6)-(III-7) that are exact at order 0{1/N). However, it seems 
difficult to go any further without either (i) considering homogeneous systems or (ii) neglecting 
collective effects. In fact, due to the huge timescale separation between the dynamical time t D 
and the relaxation time > Ntp, it could be of interest to develop a kinetic theory in angle- 
action variables as considered in [53] , where an orbit-averaged- Fokker-Planck equation has been 
derived for one-dimensional systems with weak long-range interactions (note that the formalism 
developed in the present paper could be used to better justify the kinetic equation obtained in 
|53j). This could be a future direction of investigation. Another direction of research would be 
to investigate in deeper detail the non-markovian kinetic equations derived in this paper. This 
will be considered in future works. 



A First and second moments of the velocity increment 

In this Appendix, we calculate the first and second moments (Ay* 1 ) and (Av^Av") of the 
increment of velocity of the test particle directly from the Hamiltonian equations of motion 
(1-1). We follow a procedure similar to that used by Valageas [51] in a different context. Since 
the calculations are similar, we shall only give the main steps of the derivation. For simplicity, 
we assume that all the particles have the same mass m. In order to separate the mean field 
dynamics from the discrete effects giving rise to the diffusion and to the friction of the test 
particle, we write the Hamiltonian (1-1) as 

H = m{H Q + H I ), (120) 

where we defined the mean field Hamiltonian Hq by 

H = lY. v i+Yl®°( r ^ ( 121 ) 

i i 

and the interaction Hamiltonian Hi by 



H 1 = e' 



iut 



m 

. i<3 



'1221 



In Eqs. (1 1 2 1 1) - ( 07221) the mean field potential is given by 

$o(r) = J p(i>(|r-r'|)dr', (123) 

where p(r') is the mean field equilibrium spatial density of the particles. The factor e ut has been 
added for the computation of perturbative eigenmodes and we shall ultimately let uj — > + . Thus 
H describes the mean field dynamics whereas Hi describes the discrete effects which vanish 
in the limit N — > +oo. Therefore, we consider Hj as a perturbation of H and we apply a 
perturbative analysis in powers of 1/N. We assume that the system is spatially homogeneous 
so that $o = 0. The interaction Hamiltonian Hi can be rewritten 

H I = e U)t mY^ [ w(k)e ik(r <- r ^k. (124) 

i<3 
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On the other hand, the equations of motion read 

dvi 1 dH d + H ) ^ ^ ^ (125) 
dt m dYi dri ' dt m <9vj 

We write the trajectories {r(t), v(t)} as the perturbative expansions r = r^ ^ + + r^ 2 ^ + ... 
where r^ fc ^ is formally of order k over Hi. At zeroth-order, we simply have 



which yields the rectilinear orbits 

vf(t') = Cte = v,, rfV) = Vi (t' -t) + r,, (127) 

where, in the following, r-j and v, denote the position and the velocity of the particle i at time 
t. At first order, we obtain 

*"? _ dH, dr«_ (1) 

^r- Vj ' (128) 

where we can substitute the zeroth-order orbits in the r.h.s. of these expressions. This yields 



(i) 



-jr e " tm H f Hky^-^dk. (129) 

dt dr t . <f J 

Using Eq. (I127p and integrating on time, we obtain 

v 4 (1) = -^- f dt'e^'mj^ /«(A;)e ik ^-V)(*'- t )e ik ^-v)dk. (130) 

<*i J -co j<f J 

Therefore, 

d , rf /" „ x e ik(r ^v) 

Substituting this expression in Eq. ( 1 1 28l) -b. we get 



e^ m y u{k) — -dk. (131) 



7 (1) , , ik(r (()) -r (0) ) 

^ = / «w ; :. m ' m . *. d32) 

dt w + «k( V y-v} ; ) 



Using Eq. (I127P and integrating on time again, we obtain 

■ ke ik( ri - rj ) 

/ 

This can be rewritten 

3 ik(rj-r,/) 



. . ^ 7 'L-p* k ( r i- r j) 



(ik. (133) 



9vi ^ + zk(v,-- Vi 



(134) 
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Therefore, at first order, we have 



with 

X 

At second order, we have 



r « - -A r (1) - A fl35) 



/ w(*0 n ;dk. (136) 

j<f J ^ + «k(v i -v i 



-"J 1/(1) ^oyN 



c?t drfdr. 

j 1 3 



The average acceleration of the test particle is 

(^ (2) ) = " (S^r^A - n/J^-t^ ) . (138) 



Using Eq. (I134p . the relations 

d 2 Hj 



- = -e^ m J2 J n(k)PA;V k(r - r ^rfk, (139) 



dr^dr 



J M(k)^re ik(r - ri) c/k, (140) 



and performing the averages in Eq. (11381) with respect to the distribution function /(v), we 
obtain after some calculations 

-(27T) | W( V ^«(fc)Vr ^ +( ^ w)2 . (141) 

We introduce the velocity increment Af' 1 = f M (t+ At) — f M (t). Noting that Eq. (1141ft represents 
the variation of the velocity increment at order 1/N, taking the limit uj — * + and using 
lim w ^o^/(^ 2 + % 2 ) = ttS(x), we find that 

^r) = / dvidk/(vi)«(*) W (A _ JL^) 5(k . w) . ( i42) 

On the other hand, the diffusion tensor (Av^Av u ) at order 1/N is equal to (Av^ Av "W). 
According to Eq. f 1 1 3 1 [) . we have 



ik(r(°)-r (0) ) 



-ie^'m V / ulk)^ j^dk. (143) 

f^i u; + «k(v(o)-v 0) ) 
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Therefore, using Eq. (11271) . we obtain 

e ?.k(r-r 3 )+ik(v-v J )Ai 

/ dku{k)& 

3^0 



oo + ik(v — Vj] 



r ik(r- rj ) 

+^mV / dku(k)k» —, r. (144) 



Taking the average with respect to the distribution function /(v), we obtain after some calcu- 
lations 



(2n) d me 2ujt J dkdvj^k^k 1 - 



At I v ' J 1JV iy w 2 + (k-w) 2 

x — [1 + e 2uM - 2e" A * cos(k • wAt)] . (145) 



'A* 

The limit u; — > + now gives 

{—&-) = 2 ^) dm J ^/(vOfeT [1 - cos(k • wAt)] . (146) 

Finally, taking At — > +oo and using lim^ +00 (l — costx)/tx 2 = tt5(x), we obtain 

2Ti(2n) d m / dkdv 1 f(v 1 )k fl k u u(k) 2 5(k- w). (147) 



At 

This relation can also be obtained from the Kubo formula J + °° v^^fyv^it + At)d(At) using 
Eqs. (I127p and (11291) . Equations (I142p and (I147P return the terms of diffusion and friction 
obtained in Sec. [3] when collective effects are neglected. 

B Collective effects 

When collective effects are taken into account, the force created by a particle on the others is 
modified by the influence of a "polarization cloud" . This effect can be calculated precisely in the 
case of a spatially homogeneous medium. In that case, the linearized Klimontovich equations 
are 

"'■v^-V» *=0, (148) 



dt dr dv 



<J$(r, t) = J u{r - r') [Sp(r', t) + mS(r' - n - Vit)] dr', 



(149) 



where (ri,Vi) represent the position and the velocity of the test particle (here denoted 1) at 
time t = 0, and we have made the linear trajectory approximation ri(t) = rj + V\t that is valid 
at leading order. Taking the Laplace- Fourier transforms of Eqs. (11481) and (11491) . we obtain 

5f(k, v, oo) = ^5$(k, oo), (150) 

k • v — oo 



5®{k, oo) = (2rr) d u(k) j Sf(k, v, oo) dv + m u(k)e~ i]iri 5(k ■v 1 -u). (151) 
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Substituting Eq. (IToTIj) in Eq. ffTBTD . we find that 



= m 



u(k) 
e(k,w) 



S(k • vi - ct>), 



(152) 



where e(k,u) is the dielectric function (11-50). Taking the inverse Laplace- Fourier transform of 
Eq. (11521) . we obtain the effective field created by particle 1 on particle taking into account 
collective effects: 



Finally, the effective force created by particle 1 on particle taking into account collective 
effects is 
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